Effects of Disorder in Location and Size of Fence Barriers on 
Molecular motion in Cell Membranes 
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The effect of disorder in the energetic heights and in the physical locations of fence barriers en- 
countered by transmembrane molecules such as proteins and lipids in their motion in cell membranes 
is studied theoretically. The investigation takes as its starting point a recent analysis of a periodic 
system with constant distances between barriers and constant values of barrier heights, and employs 
effective medium theory to treat the disorder. The calculations make possible, in principle, the ex- 
traction of confinement parameters such as mean compartment sizes and mean intercompartmental 
transition rates from experimentally reported published observations. The analysis should be help- 
ful both as an unusual application of effective medium theory and as an investigation of observed 
molecular movements in cell membranes. 

PACS numbers: 05.60Cd, 05.40Fb, 87.15kt 



I. INTRODUCTION 

The biophysics of cell membranes is an active field of 
current research, issues of interest being cell shaping and 
movement [1], cell division j2], signal transduction [3], 
and molecule trafficking |4]. Observations of the lateral 
movement of molecules on the surface of the cell |5J 13 
[SI [HI Uni IH] have given rise to the idea that the moving 
(transmembrane) molecules are confined within certain 
regions of the cell membrane. One possible source of this 
confinement has been suggested [TT] as being collisions of 
membrane molecules protruding into the cytoplasm with 
the cytoskeleton [1^. The model views the molecules 
as moving freely, their motion being hampered as they 
traverse adjacent compartments. As the actin filament 
that forms the compartment boundary dissociates due to 
thermal fluctuations, the moving molecules is envisaged 
as overcoming the barrier potential and hopping to the 
adjacent compartment. 

To our knowledge, two theoretical attempts have been 
made in the literature to analyze such barrier-hindered 
motion of transmembrane molecules. In both of them 
the molecule is looked upon as a random walker mov- 
ing (in a 1-D system for simplicity) with periodically ar- 
ranged semipermeable barriers. The more recent of the 
two attempts borrows the spirit of the earlier one [T3] 
but avoids some of its shortcomings. The shortcomings 
include the appearance of an unfortunately unwieldy in- 
finite series of terms which is difficult to handle, and the 
unavailability of explicit usable expressions for the mean 
square displacement, which is the quantity of direct com- 
parison to the experiment llj. It is that second (more 
recent) theoretical attempt yL4j that we take as our start- 
ing point here and calculate by substantial modifications 
of that analysis the consequences of disorder on the ef- 
fective diffusion constant and the molecular mean square 
displacement. 

In the model analyzed in ref. [14], we represent the 
molecule as a random walker in a 1-D infinite chain of 



sites. The molecule, whose probability of occupation of 
the TO-th site of the chain at time t is Pm(t), hops via 
nearest neighbor transfer rates. The rate is F within a 
compartment and has a lower value / at the interface 
of compartments where there is a barrier hindering the 
molecular motion. There are H + \ sites, equivalently 
H nearest neighbor bonds, within a compartment; for 
simplicity, H is taken to be even, with the site at the 
center of one of the compartments. Specifically, the Mas- 
ter equation 



dPrn 

dt 



F [Pm+l + Pm-l — 2Pm] 



<P - /) E [^-^+1 - P-^ ('^™.'- - '^'».-+l) 



(1) 



describes that ordered system. The primed summation 
goes over sites r = Hjl -I- {H -I- \)l which lie to the left 
of each barrier, / taking all integer values. 

Explicit expressions have been provided in ref. |14j 
for the time dependence of the mean square displace- 
ment of the molecule and for the effective diffusion con- 
stant, the compartment size H and the transfer rates / 
and F being reffected transparently in these calculated 
quantities. These are based on Eqs. (p| to appear in 
Section 2 below. The major element missing from that 
periodic barrier theory is the realistic effect of the com- 
partment sizes and the barrier heights being not equal 
throughout the system, in other words of H and / be- 
ing variable, i.e., disordered, quantities. Our interest in 
the present paper is to treat this important disorder ef- 
fect. The tool we employ is effective medium theory: a 
disordered system is replaced by an ordered system prop- 
erly structured to represent the original system. There 
is an extensive literature on the general subject of ef- 
fective medium theory [TBI HI EZl EEl EH], and there is 
a long history of its application to problems involving 
normal and anomalous transport in disordered lattices. 
Recently, the approach has been applied with success to 
the study of transport in more complex environments. 



including random graphs f^IP and small world networks 
[2T| [22l [23] , Some of those developments are naturally 
applicable to the present problem of molecular motion in 
cell membranes as we shall see below. 

The motivation to extend our previous theory [14] to 
treat disordered / and H stems from the experimentally 
known fact that the barrier locations for the transmem- 
brane molecular motion in the cellular membrane occur 
at positions that are by no means regular; the conse- 
quent variations are substantial within a system, about 
one order of magnitude, the compartment sizes some- 
times being quoted as lying between 30 nm and 240 nm 
[llj . The precise situation in which the moving molecules 
find themselves at the barriers also varies, the result be- 
ing a variation in the effective transfer rate. 

We construct our effective medium theory considera- 
tions for the present problem in three parts. First, in 
Section 2, we take the compartment sizes to be all equal 
as in ref. [13] or [T3] but allow the barrier heights, conse- 
quently the inter-compartment rates of molecular motion 
/, to be taken from each of several specific distribution 
functions with given mean, variance, and nature. In Sec- 
tion 3, we first take the inter-compartment rates to be 
constant throughout, but allow the compartment sizes 
H to vary, and then allow both H and / to be random 
variables. Fig. 1 illustrates the three respective cases as 
(a), (b) and (c). 



(a) 



Jl 



developed in ref. [U]. Fully arbitrary distributions for 
compartment size will be studied in a future publication. 
A graphical comparison of the predictions of our effective 
medium theory with the results of a numerical solution 
of the disordered problem, and concluding remarks are 
presented in Section 4. A brief analysis of the memory 
function that emerges from the effective medium theory 
is also given in the discussion in Section 2. 



II. DISORDER IN BARRIER HEIGHTS 

The system we investigate in the present paper is for- 
mally described by Eq. (IT]) as in ref. [14] but with the un- 
derstanding that / and the primed locations r are disor- 
dered quantities. In the analysis of ref. [Tlj, the solution 
of the ordered Eq. (II]) for arbitrary initial probabilities 
P„(0) is expressed as 



P.n{t) = Y.XmAt)Pn{^) 



where the transport propagator (which 'propagates' the 
solution from site n to site m) is given in the Laplace 
domain (e being the Laplace variable and tildes denoting 
Laplace transforms) by 



Xrn,n = * 



F-f 



l + {F-f)Jl 
I 



(2) 



*. 
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^lJI. 
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FIG. 1: Schematic illustration of the diflferent types of bar- 
rier disorder considered. In (a), the barriers are periodically- 
spaced, but have random energetic heights. In (b), the barrier 
heights are uniform, but the distances between barriers is ran- 
dom. In (c), both the barrier spacing and the barrier heights 
are independent random variables. 

In each case we calculate physical observables typified 
by the effective diffusion constant. Whereas the / dis- 
tributions we consider in Section 2 are arbitrary, the H- 
distributions we consider in Section are not arbitrary but 
determined by the specific / distributions we take to gen- 
erate them. This allows us to make direct use of the anal- 
ysis (in particular, the form of the transport propagators) 



This is a slightly rewritten, but completely equivalent, 
version of Eq. (2) of ref. [14. 

Here, the first term on the right hand side is the 
Laplace transform of the propagator ^m-n(i) of the sys- 
tem without barriers {F — /), and is characterized by 
a single index as a result of complete translational in- 
variance at the site level. The second term describes the 
effect of the barriers between compartments. It is pro- 
portional to the difference F — f and is characterized by 
a property of the barrier less system, viz., products of the 
propagator differences in the Laplace domain, summed 
over barrier locations, and is also characterized by fi{t), 
an appropriate summed combination of ^(t)'s (see ref. 
|14j). whose Laplace transform is given by 



^(^) = F 



tanh (^/2) 



tanh {^{H+1) /2) 



with ^ == 2sinh^ {y^e/4F). In obtaining this result, one 
uses the known form of the Laplace transforms of the 
propagators for the system without barriers. 



*/ = 



2Fsinh^' 



(3) 



A. General Development 



in the Laplace domain, 



Our goal in this Section is to analyze the generaliza- 
tion of the system represented by Eq. (pi when the in- 
tercompartmental transition rates / vary in magnitude 
throughout the chain and are picked from a distribution 
function p{f). The equation obeyed by the probabilities 
of occupation is 



F-T 



dPr,. 
dt 



— F [Prn+l + P-m-l — "^Prn] 



J2iP - fr) [Pr+1 ~ Pr] {Srn,r " <5™,r+l) (4) 



An exact analytic solution of Eq. Q is practically im- 
possible for large systems because the /r's vary in a disor- 
dered fashion. Therefore, in the spirit of effective medium 
theory [15j [TBI 113 UHl US] , we replace the actual disor- 
dered system with its many /'s by an effective medium 
system characterized by a single quantity (memory func- 
tion) ^(i) which is time-dependent and to be determined 
from the distribution p{f). The effective medium system 
is identical to the periodic system represented by Eq. 
([2J) except that / is replaced by J^{e). This replacement 
means that the occupation probabilities in the effective 
medium system obey 



dPUt) 
dt 



= F [P„+i(i) + P,n-l{t) - 2P„,{t)] 



- f dt' [F6{t - t') - T{t - t')] 
Jo 

I 
/ ^ [Pr+l{t ) — Pr{t )] {S„i^r — ^m,r+l) 



(5) 



Equation ([5| differs from the ordered counterpart (fl| in 
just one aspect: the appearance of the effective memory 
J-"(t) in place of the constant /. This memory, and the 
consequent convolution form for the equation, is an es- 
sential and well-known feature of the effective medium 
treatment. For long-time considerations, one replaces 
J-'(t) by its Markoffian approximation S{t) J^ !F{t')dt' . 
Then the formal identity to Eq. (Ill) is exact. 

To determine J- , we follow the effective medium pre- 
scription [T5J [TBJ [TTJ [T5J [TH] of considering one defect 
in the otherwise periodic system ([5]) formed by replac- 
ing T hy an f drawn from its probability distribution 
p{f), solving the defect problem exactly in the Laplace 
domain, averaging the solution over the /'s in the dis- 
tribution, i.e., carrying out an ensemble average of the 
solutions, and then requiring that the ensemble-averaged 
solution is equal to the solution of the system without 
the defect. 

The propagator for the effective medium system (l5| is. 



1 + (P - T)p 
x^($r_„+i -$^_„)($™_^-$™_^_i) (6) 

r 

which is precisely Eq. (pi with the replacement of / by 
T . For the defective system made by introducing the rate 
/ drawn from its probability distribution and placing it 
between the sites s and s -f 1, the propagator is the sum 
of the propagator given above and an additional term so 
that the defective propagator is 



l + (/-^)/? 



\X,ra,s X?n,s + 1/ vXs + l.n \s,n) 



where, for notational convenience we have introduced the 
abbreviation 

/5 = — Xs+l,s + Xs+l,s+l + Xs,s — Xs,s+1- 

The second term in the defective propagator describes 
the modification by the barrier lying between the sites s 
and s + 1 . We get a different solution for every ensem- 
ble member, the difference being in the value of /. We 
require the self-consistency condition that the ensemble 
average over the /'s give us simply Xm.n- This must be 
true whatever the n in the propagator or whichever bar- 
rier s characterizes. Therefore, the ensemble average of 
the factor in the square brackets in the propagator ex- 
pression above must vanish. This provides a prescription 
for obtaining the effective quantity T through the solu- 
tion of the implicit equation 



d!p{!) 



f-P 



1 + (/ - ^)/3 



0. 



(7) 



The chain details are reflected in /3 and the randomness 
of the /'s in p. There is no /-dependence in JF and (3, 
although /3 is a function of e as well as of J^{e). Because 
of the independent dependence of /3 on e, the solution of 
Eq. (|7| yields an explicit e-dependence of the effective 
quantity J-"(e) that we seek. Different probability dis- 
tributions result in different expressions for the effective 
medium quantity J-. 

For long times, the Markoffian approximation of the 
memory function is appropriate. This involves taking the 
e ^ limit of the Laplace transforms in the expressions 
above. We first see from Eq. ^ that, in this limit. 



P — —Xs + l,s + Xs + l,s-|-l + Xs,s — Xs,s+1 



1 



T{0) 



This remarkable simplification allows us to conclude 
that, in the long time approximation, the effective value 



of the intercompartmental rate is given trivially as the 
reciprocal of the ensemble average of the reciprocals of 
individual intercompartmental rates: 



1 



1 



feff ^(0) 



df 



f ' 



(8) 



Application of effective medium theory has thus re- 
duced the disordered problem of interest in the present 
paper into the ordered effective problem which was com- 
pletely analyzed in ref. [14] . Combining that analysis 
with Eq. ( 35 1 , we find that the overall effective transfer 



rate for the diffusion of the molecule (taking into account 
both the existence of the compartments of size H and the 
existence of disorder in the rates /) is given by 



"eff 



(^) 



ilJdf 



Pill 
f 



(9) 



For large H (for instance if _ff >> 1), Eq. ^ simply 
states that the effective overall transfer rate is the har- 
monic mean of the intracompartment rate F and the ef- 
fective intercompartmental rate feff reduced by the size 
of the compartments. 



B. Specific Cases 

Equation ([9| is one of our main results for the case in 
which disorder appears only in the values of the inter- 
compartmental rates: it allows us to translate the ran- 
domness of the intercompartmental rates as expressed in 
the form of the distribution /?(/) directly into the effective 
diffusion parameters of the system. We will now consider 
several different cases of p(/) for purposes of illustration. 

Two trivially expected results emerge in a straight- 
forward fashion: If p(/) has a non-zero value at / = 
so that there is a non-zero probability of having discon- 
nected sites, the effective hopping rate at long times must 
vanish. This is clear from Eq. (|9| since p{f) can be writ- 
ten as p{f) = 6{f) + R{f) where R{f) is some distribu- 
tion of / obeyed for all / except for / = 0. The integral 
in the denominator of Eq. [9] diverges, giving -Fe// = 0. 
Similarly, if there is only a single value of the intercom- 
partmental rate /, viz. g: p{f) = S{f — g), Eq. ^ 
reduces to the corresponding expression in the analysis 
of the ordered system treated in ref. [14j . 

For the case when there are two values of the inter- 
compartmental rate appearing with different weights: 

p{f) = A,S{f-h) + A2Sif-h), 

where obviously the normalization is Ai + A2 — 1, the 
overall effective rate is given by 



Peff - 



H+l 



/1/2 



(10) 



Note that, if either /i or /2 is zero, -Fe// vanishes as we 
have stated above. 



Finally, we will exhibit the continuum limit of our re- 
sults, selecting three specific distributions for the inter- 
compartmental rates. We display the continuum limit ex- 
pressions because, particularly for the problem of molec- 
ular motion in cell membranes, they are more directly 
applicable than their more general discrete counterparts. 
The continuum limit means that the lattice constant 
a — > 0, transforming hops among discrete sites on a chain 
to fiow on a continuous line. As is well-known, and explic- 
itly commented on elsewhere, (see, e.g., refs. [HI [24]). 
as a ^ 0, it is necessary that f,H -^ 00 as 1/a but 
F — > 00 as 1/a^. With this appropriate limiting behav- 
ior, the overall effective diffusion constant is given as 



limFe//a 

a — •^U 



D. 



D 



eff 



^/ 



1+f 

lim fa 



Jd&f 



PJS'f) 



(11) 



where D is the continuum limit of Fa^, S>f is the contin- 
uum limit of fa, and L, the continuum limit of {H+ l)a, 
is the size of the compartment. Factors such as {H+l)/H 
collapse into 1. The ratio £P — S^f/D is what is some- 
times called [13] the permeability. 

To illustrate the effect of the form of the distribution 
functions, we now consider several explicit realizations 
of p{^f). We evaluate Eq. (Ill for the three respective 
cases of a constant distribution in an interval, a biased 
distribution that peaks at a value related to the spread 
of the distribution, and a biased distribution that peaks 
at a value independent of the distribution spread. We 
use the normalization / p {^f) d^f = 1. 



1. Uniform distribution 

If p i^f) is a non-zero constant in an interval of values 
of &f, i.e., for I < ^f< u, and vanishes otherwise, then 
we have from Eq. ([Tl| 



D, 



eff 



D 



1 + 



i:>ln(u//) 
L{u-l) 



(12) 



2. Rayleigh distribution 

If p{!^f) is a biased Gaussian, called sometimes a 
Rayleigh distribution: 



P(%) 



,e-®|/2.^ 



(13) 



where the mean is ay^n/2, and the variance a^ {^^^) is 
proportional to the square of the mean, we have 



-D, 



eff 



D 



(14) 



3. Rice distribution 

To have two independently controllable parameters, 
one deciding the value at which the distribution peaks, 
and the other the spread, we consider what is called the 
Rice distribution: 



P(%) 



^/. 



a^ 



'0 



i^6) 



(15) 



the two parameters being a and v. The mean 
is ay^T:/2Li/2 (— u^/2cr^) , and the variance is 
2cr2 + 1,2 _ 7rcrV2i2^2 (-^VZcr^) . Here, ^1/2(2;) = 
e'^/^[{l-x)Io{-x/2)-xIi{-x/2)] is the Laguerre 
Polynomial of fractional order and Im{x) are modified 
Bessel Functions of the first kind. We find 



D. 



eff 



D 



-""/^''"l„{v'/4a^) 



(16) 



The distributions are plotted in Fig. 2. Note that the 
Rice distribution appears highly symmetric around its 
peak value although, like the Rayleigh distribution, it 
has the value at the value !^f — 0. 

Notice that, in every case, the expression for the ef- 
fective diffusion constant depends in essentially the same 
manner on the ratio of the system (barrier- less) diffusion 
constant D to the product of the compartment length 
and a characteristic ^/ value. Except for numerical fac- 
tors, the latter product is {u~l)/ln{u/l) for the uniform 
distribution, cr/-\/ 7r/2 for the Rayleigh distribution, and 

a/ \ ^/TTj2e-''^ /^"^ Io{v^ /Aa"^)] for the Rice distribution. 



dom been carried out in the literature. The few excep- 
tions are in the context of percolative systems in which 
the long time diffusion constant vanishes at the percola- 
tion point [16j . of the power law tail analysis given by one 
of the present authors [18^ , and prescriptions provided by 
another of the present authors [19] for the determination 
of memory functions for stress distribution in granular 
compacts. In this subsection we briefly show how to ex- 
plicitly calculate the memory function J- in the Laplace 
domain following a prescription similar to that given in 
the last of the above references [TO]. Equations (9) and 
(10) of ref. [12] should be compared to (|7| above and ( 35 1 
below in the present analysis. We recall that the prop- 
agators ^i{t) for the original problem (nearest neighbor 
rates F and no barriers) are Ii{2Ft)e~^^* and therefore 
their Laplace transforms can be written in terms of hy- 
perbolic functions as given in Eq. ([3|. This allows us to 
write 



1.15 



WAM/WHM 



r 



iK 



its 



10 



,IF 




FIG. 3: Evaluation of the memory function produced by our 
effective medium theory. Plotted is the e— dependence of the 
normalized T , taking e real for simplicity in display and nor- 
malizing it to F . Shown is the case when the distribution of 
intercompartmental rates / is a sum of two weighted delta- 
functions: the rates are either /i or /2. Our evaluation shows 
that T equals the (weighted) arithmetic mean (WAM) of the 
two rates for large e and their (weighted) harmonic mean 
(WHM) for small e, the latter representing the effective long 
time intercompartmental rate, as expected. The parameter 
values chosen for this plot are: /i — O.IF, /2 — 0.2F, H = 10, 
a = 0.5. 



FIG. 2: Plots of Rayleigh (solid) and Rice(dashed) distribu- 
tions for various parameter values. For the Rice distribution, 
(T=0.05 in all cases. Distributions are normalized such that 
Jg°° dxp{x) = 1. Here, x represents ^/. 



C. Determination of the Memory Function 

The full exploitation of the convolution in Eq. (|5l and 
of the consequent memory effects in the motion has sel- 



9-^ 



F^{l + {g-T)Jj) 



/3 = 2(*i-*o) 



1 - 2e^o + e 



E^ 



2 

s—r 




cothC(iy-M) 
4F2 sinh^ $ 



(17) 

(18) 
(19) 



In the light of this expression for /3, Eq. n7[ yields: 



A. Disorder in Barrier Placement Only 



where 



r J ' f + iv-T) 



(1 + JIF) - TJl 
C + 0((1 + JIF) ~ TJl) ' 



(20) 



(21) 



In order to write the expressions in a compact way, we 
have defined: 



coth(^/2) - 1 
F 



c, 



C = :=^^(1 - 2coth(e/2) + coth2(C/2) coth^(i/ + 1)). 
ImF'' 

(22) 

We can now solve for T for a given p(/) by using Eq. 
( 20 1 as outlined for a different case in ref . [12] • As a 
special case of the distribution we use the case when the 
intercompartniental rate takes on one of two values with 
different weights: 

pif) = aSif - h) + (1 - a)5{f - /2). (23) 



Then Eq. (20) gives us, for T, 



T^ + bT^- + cT + d= 0, 



(24) 



where 



h = -(/i + h- i/e) - v/oji, 

C = hh - [m(/2 + a(/l - /2)) + (1 + llF) 

- vifi + f2)]/eji, 

d = N/1/2 - (/2 + a(/i - /2))(1 + liF)]/9Jl, 
rj^iC + eil + liF)). (25) 

The numerical solution of the cubic equation yields the 
memory function explicitly in the Laplace domain. It is 
plotted in Fig. 3. We see that, at e = 0, J^ tends to 
the value of its Markoffian approximation which is the 
weighted harmonic mean of the two rates /i and /2 . We 
also see that it tends to the weighted arithmetic mean 
of the two rates as e —^ 00. The intermediate behavior 
corresponds to intermediate times. We defer to a forth- 
coming publication [25] a detailed application of these 
and related memory developments. 



III. DISORDER IN BARRIER PLACEMENT 
AND IN INTERCOMPARTMENTAL RATES 

We now turn our attention to a generalization of the 
problem which involves disorder in the placement of the 
barriers rather than in their heights. We first take the 
latter to be the same throughout the chain. 



Let us consider the following distribution for barrier 
heights: 



p{f)^adif-g) + il~a)6if-F). 



(26) 



Here, the intercompartniental transition rate is either g 
or F with probabilities a and 1 — a respectively. When 
the barrier heights are distributed according to Eq. ( 26 1 , 



it is as if, starting with exactly the same system as in ref. 
[14] (periodic barriers of the same height g) , we replace 
the barriers randomly and independently with probabil- 
ity 1 — a with links that have transfer rates F, which 
is the intracompartmental transfer rate. Therefore, the 
barriers between some of the compartments are now re- 
moved which means that some of the compartments are 
merged. The compartment sizes always are in multiples 
of {H + l)a (the smallest compartment size) where a is 
the lattice constant. In this system, the compartment 
size, which is the distance between the consecutive bar- 
riers, will be a random variable whose distribution will 
depend on the parameter a. 

The distribution of distance between two consecutive 
barriers is easy to calculate. Consider TV points in a dis- 
creet linear space. To each point we will assign a number 
Si which is either 1 or 0. We will think of the points in 
this space as the intercompartniental links in our original 
problem. In this picture, a point i with Si — 1 will repre- 
sent a link with transfer rate F and a point with Si = 
will stand for a link with transfer rate g, which we have 
been calling a barrier. Then according to Eq. ( [26| ), O's 
will occur with probability a and I's with 1 — a. If we let 
a be the number of elements in a contiguous sequence of 
I's, the distance between two consecutive barriers is sim- 
ply given by ((7+l)(iJ+l)a. Note that ct = corresponds 
to the case in which the distance between two consecu- 
tive barriers is the maximum value (H + l)a. When all 
these arguments are taken into account, the number dis- 
tribution of a is found to be given by 



N-l 
J = l \i=0 ) 

(27) 

where y^{(7) is the number distribution of cr in a par- 
ticular realization of a 1-D chain as described in the be- 
ginning of this section. The first and second terms in 
Eq. ( 27 1 count the occurrence of compartments of sizes 



{H H- l)a and {a + 1){H + l)a respectively. As s^'s are 
independently distributed, we can write: 



{S^) 



1 — a, 



(28) 




FIG. 4: The ensemble averaged probability distribution 
P{a,a), in the limit as A'' — > cxa, as a function of a (left) 
and as a function of its parameter a (right). The continuous 
lines in the plots correspond to our formula in Eq. (|32|). 



above, with a prescription to calculate the effective long 
time transfer rate F^ff in terms of the mean (dimension- 
less) compartment size q and the rates F and /: 



Peff - 



l/f+{q-l)/F- 



(35) 



In the continuum limit, obtained by multiplying i^e// 
by a^ and letting a tend to zero appropriately, we note 
that H also tends to infinity, producing the limit of qa 
as the mean compartment size Q which has dimensions 
of length. We get 



D 



eff 



D 



1 



D 



(36) 



where the angular brackets mean an ensemble average 
over all realizations of 1-D chains with intercompartmcn- 
tal transition rates sampled from Eq. (26 1. Then the 



ensemble averaged number distribution is given by 

i^Ji{G)) = {N -a-\)c?{\~a)\ (29) 

and therefore the probability distribution for a is 



(^(^)> 



El^o'(-^(^))' 



(30) 



As we are interested in particularly in infinite chains, we 
take the limit N —f oo. The probability distribution for 
a becomes 



(Pw^ooM) = «(!-«)'. 



(31) 



Then the ensemble averaged compartment size distribu- 
tion would be given by 



P{a,a) = a(l - a) "+1 \ 



(32) 



where q — {a + l)(i? -I- 1) is the dimensionless compart- 
ment size. The mean and variance of P((T, a) are given 

by 



H+l 

a 



(33) 



{Aqr^{q^)^{qf^iH + ir^. (34) 



In Fig. |4]we display this distribution: plotted against 
a for three values of a on the left side and against a for 
three values of a, as shown. The left plot shows that, 
in all cases the distribution is peaked at vanishing com- 
partment size. The right plot shows that, with the ex- 
ception of the case a — 0, the distribution vanishes for 
both extremes a — and a = I, but rises and drops for 
intermediate values. 

The general theory we have developed in Section II now 
provides us, in light of the distribution we have obtained 



B. Disorder also in the Intercompartmental Rates 

We now give expressions for the effective hopping rate 
and diffusion constant when both the heights and places 
of the barriers are random. 

Consider 



p{f) = {l-a)S{f-F) + Tj{f,a), 



(37) 



where ri{f, a) is a distribution normalized to a, with the 
understanding that 77(0,0;) — 0. According to Eq. (37 1, 
and the development in Section |III A| a fraction a of 
the intercompartmental links are barriers whose heights 
are sampled from the distribution ?/(/, a) and the rest 
are just intracompartmental links with transition rates F 
that in turn give rise to the variability in compartment 
sizes. Note that the statistics of different compartment 
size distributions do not change even if the barrier heights 
are not the same. Therefore the compartment size dis- 
tribution can still be obtained from Eq. (31). Thus we 
get 
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(38) 



and the effective diffusion constant in the continuum limit 
becomes 
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(39) 



Here, /„ dxrj^x, a) = a, x being / and &f in the two 
respective equations above. Note that when a = 0, so 
that there are no barriers, ri(x, 0) vanishes identically as it 
is a positive function, and the results reduce to F^ff — F 
and Deff = D. 



IV. CONCLUSION 

This paper has a twin purpose, an investigation of 
the effect of disorder on molecular motion in cell mem- 



branes, and the development of effective medium the- 
ory approaches in new practical directions. We have dis- 
cussed the first context briefly in the Introduction. Sin- 
gle fluorescent video imaging [2SJ H7] and single parti- 
cle tracking [28i 2^ are the two most common ways of 
observing laterally moving molecules on live cell mem- 
branes. The latter can provide information about mo- 
tion at very short times such as confinement effects. Our 
theory is positioned to address the results of both mea- 
surement techniques. Our results allow us to describe 
the effect of disorder both in the location and the size 
of fence barriers, and the motion of both proteins and 
lipids in cell membranes. These barrier parameter vari- 
ations arise from the dynamics of the cytoskeleton and 
of the transmembrane molecules. Observational deter- 
mination of such variations should be possible through 
the use of modern methods involving optical tweezers 
[30] and electron tomography techniques Zl\ that enable 
the imaging of the actin filaments and the cell surface. 
These experimental data can be used in conjunction with 
the calculations we have presented to extract informa- 
tion about compartment sizes and their spatial fluctua- 
tions, as well as about the energetic barrier magnitudes 
and their fluctuations. As more detail becomes available 
about the membrane dynamics from experiment, it will 
be straightforward to give further quantitative shape to 
the theory we have provided here. 

The second context of our study lies in presenting ef- 
fective medium calculations in systems which are partly 
ordered and partly disordered. This line of research has 
been recently taken by two of the present authors in their 
study of transport on small world networks, particularly 
of the Neumann- Watts kind [21j i22i i23j . In those systems 
standard rings (finite chains with periodic boundary con- 
ditions) with nearest neighbor hopping rates for the ran- 
dom walker form the ordered part and additional small 
world connections make up the disordered part . Effective 
medium theories developed for those systems envisage 
the effective system as the rings with additional period- 
ically placed connections. Here, for the cell membrane 
problem, the standard chain forms the ordered part, and 
the disordered barriers (in magnitude and location) the 
disordered part. Our effective medium theory here takes 
the effective system as the chain with periodically placed 
barriers of constant magnitude. 

Our starting point in this paper has been Eq. Q which 
is impossible to solve in practice because of the enor- 
mous number of irregularities. Effective medium consid- 
erations have yielded as a general result the prescription 
( 7 1 to calculate T appearing in the effective medium Eq. 
( 5 1 as a function of the Laplace variable e, starting from 
known random distributions. We have illustrated the 
procedure to carry out this prescription for a si mple case 
which results in a cubic equation for T, viz. (24 1. Its 
solution is in Fig. l3] and shows physically expected lim- 
iting results. Long time approximations to the memory 
!F{t) are appropriate if the interest lies in long time de- 
scription of the transport of the random walker (e.g., of 
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FIG. 5: Excellent agreement for sufficiently long times of our 
effective medium theory with numerically obtained solutions 
of the Master equation with disordered barrier heights typ- 
ified by three Rice distributions with different parameters, 
{v = 0.1, s = 0.02), {v = 0.2, s = 0.06) and (v = 0.4, 
s = 0.02), respectively denoted by dashed, solid and dash- 
dotted curves in both the main figure which depicts the in- 
stantaneous transfer rate F{t) normalized to F, and in the 
inset which shows the corresponding probability distributions 
p. Long time predictions (effective transfer rates) of our ef- 
fective medium theory, shown by horizontal dotted lines, are 
reached asymptotically in each case. See text for other details 
including the observed dips. 



the molecules in the cell membranes.) In such a case we 
have developed a quite general result, Eq. (Ill, which 



transparently connects the distribution of intercompart- 
mental rates p{f) to the effective diffusion constant. We 
have illustrated these results for several specific cases of 
the distribution p(/) including a Rayleigh and a Rice dis- 
tribution, and also employed the general development to 
shed light on effects of disorder in barrier placement as 
well as in intercompartmental rates. 

In concluding, we display a comparison of the results 
of our effective medium theory with numerical solutions 
of the disordered Master equation for various realizations 
of /-disorder, with H held constant. We see in Fig. [5] 
that the agreement is excellent at sufficiently long times. 
We considered several different kinds of distribution and 
found the results to be essentially identical. The results 
displayed in Fig. [Slcorrespond to three Rice distributions 
with different parameters, (u = 0.1, s — 0.02), {v = 0.2, 
s = 0.06) and {v = 0.4, s = 0.02), the probabihty dis- 
tribution functions being shown in the inset. The same 
three kinds of curve, dashed, solid and dash-dotted curves 
respectively, used to display the distributions are used 
correspondingly in the main figure to show the results of 
the numerical solution of Eq. (Hi) followed by performing 
an ensemble average. The specific plots are of the instan- 



taneous transfer rate F{t), normalized to the barrier- less 
system transfer rate F. This F{t) is one half the time 
derivative of the (dimensionless) mean square displace- 
ment, and should not be confused with the memory J^{t). 
Two features are visible: dips below the eventual asymp- 
totic values, and the coincidence of the asymptotic values 
with horizontal dotted lines that represent our effective 
medium theory. The dips arise from the fact that the 
random walker is assumed to start initially at the center 
of one of the compartments: repeated encounters with 
walls when the effective transfer rate drops are responsi- 
ble for the dips. A detailed discussion of the dips is given 
in ref. |14] . The asymptotic coincidence of the numeri- 
cal solutions with the effective medium theory provides 



graphical validation of the latter. 
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